Trigonal Band Structure and Time-Reversal Invariance in Graphene 
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We present a symmetry analysis of the trigonal band structure in graphene. While the energy 
spectrum near the Fermi edge equals the spectrum of massless Dirac fermions, the transformational 
properties of the underlying basis functions are qualitatively different. Using group theory we de- 
velop an invariant expansion of the Hamiltonian for the electron states near the K points of the 
graphene Brillouin zone. We find that the fc-linear dispersion near the band edge arises as an unusual 
consequence of time-reversal invariance. We suggest to divide the electronic properties of graphene 
into two categories, those that depend and those that do not depend on the transformational prop- 
erties of the Bloch functions at K. 
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In recent years, tremendous interest has been focused 
on studying single layers of graphene, following the first 
experimental realization of this material [l|. To a large 
extent, these continuing efforts are motivated by the 
unique band structure of graphene in the vicinity of the 
Fermi edge. In most semiconductors, the band edges 
are characterized by a quadratic dispersion, with fc-linear 
corrections possible only in inversion-asymmetric mate- 
rials due to spin-orbit coupling ^. In contrast, for 
graphene, the dispersion £'(k) of the uppermost valence 
band and the lowest conduction band is dominated by k- 
linear terms j3i], E(]i.) « hvk with Fermi velocity v. These 
bands touch at the points K and K' at the edge of the 
Brillouin zone [Fig. [Ijb)] so that the resulting energy sur- 
faces resemble those of free massless fermions described 
by the Dirac equation [Fig. [ijc)]. The apparent analo- 
gies between a solid-state system and relativistic quan- 
tum mechanics have greatly stimulated the interest in 
graphene \M,\M- I^O'^ example, it was argued that the chi- 
ral electrons in graphene can show both weak localization 
and antilocalization (Ref. j5i] and references therein). 

The group-theoretical analysis presented here shows 
that, nonetheless, there are important differences be- 
tween electrons in graphene and massless Dirac fermions. 
Even though both systems can be described by the same 
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FIG. 1: (a) Honeycomb lattice of graphene. Atoms in sublat- 
tice A (B) are marked with open (closed) circles, (b) Brillouin 
zone and its two inequivalent corner points K and K'. The re- 
maining corners are related with K or K' by reciprocal lattice 
vectors, (c) Dispersion E{k) near the K point. 



Hamiltonian (giving also the same spectrum), the trans- 
formational properties of the underlying basis functions 
are qualitatively different. We suggest to divide the elec- 
tronic properties of graphene into two categories, those 
that emerge from the linear dispersion but that are in- 
dependent of the transformational properties of the basis 
functions and those that do reflect these transformational 
properties. 

The starting point of our analysis are the Bloch func- 
tions at K and K . In the commonly used tight-binding 
(TB) model of Refs. [3|, |6| for the bonding and antibond- 
ing TT bands in graphene, each of the Bloch functions can 
be chosen to be nonzero only on sublattice s = A or B 
[Fig. ma)], i.e., 
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Here (j>Tr{r) denotes the carbon tt orbital and the sum 
runs over the atomic positions R^ in sublattice s. The 
group of the wave vector K is isomorphic to the trig- 
onal point group D^h (while the point group of the 
honeycomb structure is Dqu)- Projection of ^K,s(r) 
on the irreducible representations (IRs) of D^h shows 
[2, \M, that these functions transform according to the 
two-dimensional IR Fg [8|, |9|]. More specifically, un- 
der the symmetry operations of -Da/i, the Bloch func- 
tion ^K,A(s)(r) transforms like |p_(^)), where \p±) = 
-j=\x:kiy). Thus, from a symmetry point of view, we 
may identify {^K,s(r)} with {|pip)}. Note that ^k',s — 
^K s [i-^-' the existence of two inequivalent valleys K and 
K' does not cause extra degeneracies, see Eq. ^ below]. 
The role of \p±) is thus reversed at the K' point: |p_) 
(|p+)) corresponds to sublattice B {A). 

The TB model of Refs. @, Q neglects the spin degree 
of freedom. Thus we have obtained an ordinary IR (Fg of 
D^h), which differs qualitatively from the double-group 



(spinor) IRs characterizing a particle with a genuine spin- 
1/2 degree of freedom coupled to a particle's orbital mo- 
tion. The double group Gd corresponding to a group G 
can be written as G^ = G ® -EG, where i5 is a rotation 
by 27r around an arbitrary axis. The presence of E in Gd 
reflects the well-known fact that spin- 1/2 spinors (which 
may be used [s!] as basis functions for the spinor IR Fy 
of D^h) change sign when rotated by 2tt. The basis func- 
tions of Fg do not change sign when rotated by 2tt, i.e., E 
acts like the neutral element E. Furthermore, the point 
group D^h contains three "vertical" reflection planes a^, 
as symmetry elements [8|. Unlike spin-1/2 spinors, the 
basis functions of Fg change sign under these reflections. 
To proceed, we use the theory of invariants [21 to con- 
struct an invariant expansion for the effective Hamilto- 
nian ?i(/C) describing the electron states in the vicinity of 
K and K'. Here K represents a general tensor operator, 
which can depend, e.g., on the components of the kinetic 
wave vector k and external electric and magnetic fields 
£ and B. K. can be decomposed into irreducible tensor 
operators KS*^' ' that transform according to the IRs F^ 
of D'ih. Using the coordinate system in Fig.[T]^b), we get 
the lowest-order tensor operators lO'^' ' in Table H] Fur- 
thermore, Ti{K) can be decomposed into blocks TiapiJC), 
where a and /? denote the spaces of the Ua- and un- 
fold degenerate band-edge Bloch functions, which trans- 
form according to the IRs F^ and Tp of D^h. (Without 
spin, we have only one band Fg. Below, we also derive 
the multiband Hamiltonian that takes into account spin.) 
For each block HafsiJC), one can find a complete set of 
linearly independent n^ x n^-dimensional matrices ^; 
that transform according to the IRs F^ (of dimension L^) 
contained in the product representation Fq x Fo. Choos- 
ing basis functions for Fg that transform like {\p-), \p+)} 
we obtain the basis matrices listed in Table HIl Then each 
block TiapiJC) can be written as 
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with material-specific coefficients a"(; 



_^^. By construction, 
each block TiapifC) is invariant under the symmetry op- 
erations in D^h- Exp licitly we get up to second order in k 



(for B = £ = 0) 
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1 + allik^a^ + kycjy) + a\l{kl + kl)t 



+ 42[{kl - kl)a^ + 2k^ky<jy\ + 0(fc3) . (3) 

Third-order terms can be read off from Tables U and [III 
Terms of yet higher orders are constructed in a simi- 
lar way. Equation ^ applies to K. To describe K' 
within the same coordinate system, we need to substi- 
tute (ux, Uy, Oz) — > ((Tx, ^fj/i ^fz) (throughout Table HI)) . 
7i|g gives rise to a dispersion 
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TABLE I: Irreducible tensor components for the point group 
Di,h- (Terms proportional to kz are not listed as they are 
irrelevant for graphene.) Notation: {A, _B} = \{AB + BA). 

i 1 i; rZj; H~ Ky ; \n,x •, orCy — i^x S'l ^x^x ~r rZyOy] 

t 3 J^xf^x \ i^yi^y 

i 4 JDxKy JDyfCx'f C-z 

J- 5 J^x 7 J^y 1 -Dy rCy — J~>x f^x 7 J^x r^y \ DytZx 5 K,yOz •, ' f^xOz 

i 6 kix^ K,y] \r^y \ rZxi f^y "~ f^x J ^ ■^\n^x 7 f^y Ji 

\r£xf i^x ~u rCy I, \k,y, Kj. + rvy f] JDzf^y, J^z^^x ] 
^x-, ^y ) i^y^y f^x^x^ f^x^y r rCyC^x 



TABLE IL Symmetrized matrices for the invariant expansion 
of the blocks Ha/s for the point group Dsh- 



Block 


Represent at ions 


Symmetrized matrices 


Hee 


Fe X n 


Fi 


1 




= ri +r2 + r6 


Fa 


O-z 






Fe 


(Jx, CTy 


Ti-ss 


FsxP^ 


Fi 


1 




= Fi + Fa + Fs 


Fa 


0"z 






F5 


""a:, ~Cfy 


W99 


Fg X n 


Fi 


1 




= Fi + r2 + Fa + F4 


F2 


Cz 






Fs 


o-x 






r4 


ay 


■^89 


Fg X r<:s 


Fs 


1,-icrz 




= Fs + Fe 


Fe 


crx,ay 



model |6|1. We get a^o = <^2p 
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We find explicit expressions for the prefactors a^^ 
by comparing Eq. Q with the dispersion of the TB 
66 - .^ all = ^a{t ~ se2p), 
^^^ V- "-^p;> ""- -62 - |a^(i- se2p), where a 
is the lattice constant, e2p the on-site term, s the overlap 
integral, and t the transfer integral [GJ]. 

The part of Tigg linear in k is formally equivalent to 
the Dirac Hamiltonian for massless (chiral) fermions [13| . 
Also, it is related via a simple unitary transformation 
with the Dresselhaus 1J| and the Rashba ^5,] term in 
quasi-2D systems [l6| . All these models give rise to a dis- 
persion that is linear in the limit of small k. Yet for each 
of these models, the transformational properties of the 
basis functions under the corres pon ding symmetry oper- 
ations are qualitatively different [17| . The Dirac equation 
reflects the continuous symmetries of the Lorentz group 
[l8|. In quasi-2D systems, the fc- linear Dresselhaus term 
is intimately related with the tetrahedral symmetry of 
the zinc blende structure [ij] , whereas the Rashba term 
emerges from a model with axial [16|] (or hexagonal [19| ) 
symmetry. These terms refer to electron states trans- 
forming according to spinor representations of the cor- 



responding crystallographic point groups, and they are 
nonzero only as a consequence of spin-orbit (SO) cou- 
pling [l6j. In contrast, the Hamiltonian Hqq refers to the 
basis functions (H]). It is applicable to spinless particles 
or particles for which the spin degree of freedom is decou- 
pled from the orbital motion, and the group element E 
does not play a role [20|. From a more general perspec- 
tive, we see here that the symmetries of a system deter- 
mine the invariant expansion ^ and the band structure 
-E„(k). Yet it is, in general, not possible to follow the 
opposite path and infer the symmetries of a system from 
the Hamiltonian and the band structure _E„(k). 

Time-reversal invariance results in additional con- 
straints for the allowed terms in the invariant expansion 
^. The valleys K and K' are inequivalent points in the 
star {k} that characterizes the IRs of the space group for 
these values of k. When we apply the general analysis 
in Ref. [4I to graphene, we find that the bands at K and 
K' belong to the case denoted "02" in Ref. [4I, i.e., under 
time reversal 9 the basis functions {ipi} at K and K' are 
linearly related via a unitary matrix T 



6'V'K,z = i'K,^ = Xl^J^^K'j 



(5) 



For our basis {|p=p)} we get T = era, (consistent with 
the discussion above). Equation ([5]) results in the condi- 
tion ^] 

T-'-HUR-'}C)T = niMK-) - KMK-) (6) 

for the diagonal blocks Tiaa{K-). Here, * denotes complex 
conjugation and t transposition, R^^ is an element of 
D^h that maps K on K (e.g., inversion), and / depends 
on the behavior of K. under time reversal, k and B are 
odd under time reversal so that then / = — 1, while £ 
has / = 1. Then Eq. ^ implies that all terms in Eq. 
^ are allowed by time-reversal invariance, yet, e.g., the 
third-order invariant ky(3k^ — ky) az is forbidden. 

The existence of fc-linear terms in Tigg is intimately re- 
lated to the behavior of Tigg under time reversal charac- 
terized by Eq. ([6]). For comparison, consider the F point 
k = of a material with point group D^h, where Tables H] 
and [ni are valid, too. Yet the fc-linear terms in Tigg are 
then forbidden because in this case the matrices {X^'} 
can be classified as even or odd under time reversal with 
{Xj ' } even ^j . We remark that the above analysis based 
on Eqs. (O and ([6]) is not compatible with a continuous 
SU{2)-\ike "valley isospin" degree of freedom proposed 
in previous work [5J, |2]J . 

It follows from Tables |T] and |TT] that a field £z does 
not affect the electrons' motion up to terms linear in 
k. The lowest-order invariant involving electric fields 
reads a23{kx£y — ky£x)<yz- Similarly, a linear Zeeman 
splitting in a field B|| = {Bx,By,0) (and independent of 
k) is forbidden by symmetry. Yet the allowed invariant 
Hgg = a^Bzdz implies that electrons in graphene have 



an orbital magnetic moment in z direction. In contrast, 
for neutrinos, which are almost massless Dirac fermions, 
the magnetic moment is proportional to their mass [22l |. 
The spin magnetic moment of Rashba and 2D Dressel- 
haus electrons has both a z and an in-plane component, 
though generally these are different [16|. 

In a field i3z > we can utilize the usual [16[ ladder 
operators a* — Isikx ± iky)/^/2 with magnetic length 
Ig — yJhl\eBz\ to obtain up to linear order of k and Bz 
\a± = \{a,j, ±iOy)\ 
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a^cTzp + a (T±) ±a22Bz<Tz , (7) 



where the upper (lower) sign applies to K (K ). Both 



Ti-Qg and Tigg have the same Landau spectrum 
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with positive integers n. Finding Eq 7^ afp signals broken 
particle-hole symmetry. For a22 = 0, Eq. ([8|) is identical 
to the spectrum of 2D massless Dirac fermions in a mag- 
netic field Bz [231 . The spectrum ([5]) is obtained also for 
the Rashba [15!| and 2D Dresselhaus [l6| models in the 
limit of an infinite effective electron mass. This result 
illustrates the fact that the spectrum ([8]) for Bz > and 
022 = is determined by the dispersion ^ a,t B — 
but does not depend on the transformational properties 
of the underlying basis functions (which are different for 
these models). 

The Hamiltonian ([3|) is spinless. Adding a spin-1/2 
degree of freedom {|T), |i)} to the basis functions of Fg 
gives rise to spinors that transform according to the 
double-group IRs Fg (with representative basis functions 
{\P- T), |P+ i)}) and Fg {{\P- i), \P+ T>}) of D,„. The 
corresponding basis matrices are also listed in Table [Til 
Thus we obtain the 4x4 Hamiltonian 



n 



Hs8 Ti-. 
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where each block Tiap is given by an invariant expansion 
of the form ([2]). For the diagonal blocks, Eq. ([6]) again 
applies. (The IRs Fg and Fg also belong to the case a2.) 
The most important consequence of SO coupling is the 
matrix element A — „ ^^ 2 (a^|[Vy x pjzly), where V is 
the microscopic crystal potential of graphene, so that we 
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which implies 



the opening of a gap A between the Fg and Fg states 
2J]. This gap is analogous to the gap A that separates 
the topmost valence band in semiconductors like Ge and 
GaAs from the spin split-off valence band [la]- Recent 
calculations suggest that A « 1 /zeV 25| . The sign of A 
cannot be infered from our analysis. 



Unlike the case of inversion-asymmetric crystal struc- 
tures like zinc blende and wurtzite, SO coupling in inver- 
sion symmetric graphene does not give rise to spin split- 
ting, i.e., with spin taken into account we get a two- fold 
spin degeneracy for all bands throughout the Brillouin 
zone [2| . Accordingly, the Hamiltonian ^ preserves the 
two-fold spin degeneracy (for B = £ = 0). An external 
field 6 breaks the inversion symmetry and Tables HI and HIl 
show that, to lowest order of k, the resulting spin splitting 
is due to the invariants o^^ {kxSy — ky£x)(Jz and several 
off-diagonal invariants [^ 



because of SO coupling, whereas 033 



The latter are nonzero only 
afa stay finite 

in the limit of vanishing SO coupling. Thus we expect 
the diagonal invariants due to fields £\\ to be dominant. 
Finally, we want to comment on the significance of 
the different transformational properties characterizing 
the wave functions of electrons in graphene and mass- 
less Dirac fermions. Our group-theoretical analysis sug- 
gests to divide the electronic properties of these sys- 
tems into two categories: (I) those that emerge from 
the linear dispersion (|4]) but that are independent of the 
transformational properties of the basis functions, and 
(II) those that do reflect these transformational proper- 
ties JJj . Clearly, the experimentally verified |27l] Landau 
spectrum ([5]) belongs to category (I). Similarly, Zitter- 
bewegung-\ike effects [28|, i.e., phenomena arising due to 
the interference between electron states from neighboring 
bands [29|, generally belong to category (I). In contrast, 
the magnetic moment belongs to category (II). Weak 
(anti)localization reflects the behavior of electrons mov- 
ing along time-reversed paths [30|] , which is qualitatively 
different for particles with integer and half-integer spins. 
Thus weak (anti)localization belongs also to category (II) 
and our analysis implies that within the framework of 
the spinless Hamiltonian ^ graphene shows weak local- 
ization. More generally, the formation of multi-particle 
multiplets belongs to category (II). Similarly, Kramers' 
degeneracy only holds for particle states transforming ac- 
cording to spinor IRs. In this context it is not important 
whether the symmetry group is discrete or continuous as 
in both cases we can distinguish ordinary and spinor IRs. 
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